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DELTA METHOD, AN EMPIRICAL 
DRAG BUILDUP TECHNIQUE 

Richard C. Feagln and William D. Morrison Jr. 

Lockheed-California Company 
Burbank, California 91520 


SUMMARY 


The Lockheed-California Company, under NASA Ames Contract No. NAS2-8612, 
Mod 6, has applied empirical drag correlation techniques to 19 subsonic and 
supersonic military aircraft and 15 advanced or supercritical airfoil concepts 
to develop an empirical drag estimation technique which can be applied in 
future design activities. The resulting method is presented in this report. 
The use of the technique provides a capability of estimating the total con- 
figuration drag polar near the cruise lift coefficient (Design C^±0.30) 
and a speed range from Mach 0.40 to approximately Mach 2.0. Included also 
is the capability of predicting the subsonic off-design performance of 
advanced or supercritical airfoil sections. Buffet onset may also be esti- 
mated. The method can be applied to wind tunnel models as well as to full 
scale configurations. The technique has been converted into a computer 
code which is compatible with the NASA Ames computer facilities. The 
program, "Empirical Drag Technique (EDET)", is presented in reference 2. 
Results obtained using this method to predict known aircraft characteristics 
are good and agreement can be obtained within a degree of accuracy judged 
to be sufficient for the initial processes of preliminary design. 
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INTRODUCTION 


In a conceptual and preliminary design atmosphere, it is dcisirable 
to employ quick methods of configuration evaluation in t.he first step 
towards configuration selection, and necessary when there are numerous 
configurations under consideration* The emperical approach of analysis 
is to use information already known about existing aircraft to predict 
the characteristics of future designs. Such an approach is adopted for 
this study. 

The Lockheed-California Company, under NASA Ames Contract No. NAS2-8bl2, 
Mod 6, has applied the empirical total drag technique of reference 1 on 
19 supersonic and subsonic military aircraft and 15 advanced or super- 
critical airfoil concepts to develop an empirical drag estimation technique 
which can be applied in future preliminary design activities. Correlated 
estimations of design lift coefficient and Mach number, compressibility and 
pressure drag, and buffet onset were generated, and the resulting method 
is presented in this report. 

This method can be used to estimate the total configuration drag polar 
within a C^ range of 0.0 to 0.60 or buffet onset and a 5»peed ratige of 
0.20 < M < 2.0. Included also is the capability of predicting the subsonic 
off-design performance of advanced or supercritical airfoil sections. Buffet 
onset may also be estimated. The method may be applied to wind tunnel 
models as well as to full scale configurations. 

The technique is presented in an easily followed format, and example 
calculations are presented. Results obtained using this method to predict 
known aircraft characteristics are good and agreement can be obtained within 
a degree of accuracy sufficient for establishing trends and for the initial 
processes of preliminary design. 

This drag prediction technique, has been programmed into a computer 
code for use on the NASA Ames IBM 360 and CDC 7600 computer facilities. 

The program, ^’Empirical Drag Estimation Technique (EDET)”, is presented 
in reference 2. 


To provide the basic data for this study, lift and drag data, buffet 
boundary, and geometric configuration information was supplied by NASA Ames 
on 18 subject aircraft. These aircraft were the T-2B, T-37B, KA-3B, A-4F, 
TA-4F, RA-5C, A-6A, A-7A, F-4E, F-5A, F-8C, F-llF, F-lOO, F-101, F-104G, 

F-105B, F-106A, and XB-70. Data for the S-3A and 15 advanced concepts 
was supplied by Lockheed. 

The data on the subject aircraft Included: 

• Wing geometry (aspect ratio, sweep, thickness, and airfoil 
section) 

• Cross-sectional area distribution (area progression curves) 

• Drag coefficient variation with lift coefficient and Mach number 
Cpolars) 

• Lift coefficient buffet limits versus Mach number 
Items correlated for each aircraft Include; 

• Design lift coefficient 
s Design Mach number 

• Drag divergence Mach number 

• Compressibility drag versus Mach number 

• Pressure drag versus Mach number 

• Lift coefficient for buffet o^et versus Mach number 

The drag prediction technique resulting from the above correlated items 
is discussed in Section 1 and the correlations themselves are presented in 
Section 3. Examples of the use of this technique for drag prediction and 
data sensitivity are given in Section 4. Basic data for individual aircraft 
are given by references 3 through 22 and the resulting data packages from 
which the correlations were obtained are collected in reference 33. 
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1. DRAG PREDICTION PROCEDURE 

The Delta Method, drag correlation technique can be used to predict 
drag "polars" for subsonic and supersonic aircraft configurations from the 
limited geometric parameters known during the advanced stages of design 
work. This method is uniquely capable of predicting the design and off- 
design performance of advanced airfoil sections and supercritical wings. 

It is most applicable near the cruise condition and is not appropriate for 
analysis of turning capabilities or performance near maximum lift coefficients. 
The technique can be applied equally well for total aircraft configurations or 
analysis of wind tunnel data. 

The technique is to be used in a preliminary design atmosphere where 
design details are minimal and where it is usually desirable to investigate 
the trade-off of such design parameters as sweep, aspect ratio, thickness, or 
body geometry on performance. All that is needed to provide input to this 
procedure is a three-view drawing or sketch of the basic configuration and a 
rough estimation of the proposed fuselage area distribution. With these data 
for input, a matrix on configurations can easily be analyzed and design 
trends noted in a time span considerably shorter than that required of 
most conventional methods. 

As with most empirically derived procedures, the range of applicability 
of this method is most accurate within the range of data from which it was 
derived. Where such limitations exist, they are pointed out in the following 
discussion so that caution can be exercised if the user wishes to analyze 
configurations outside this range. The computer program of reference 2 is 
also structured so that these limits cannot be passed thereby eliminating , the 
inadvertent calculation of data outside the region of applicability. 

It should be noted that the drag correlations to follow do not Include 
the effects of maneuvering devices such as flaps, slats, or boundary layer 
control and is representative of a clean wing configuration only. The effects 
of high lift devices if desired, will have to be added to the final results of 
this method by the individual concerned using conventional estimation procedures. 
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The flight test drag polars for all subject aircraft were broken down 
into their individual components by the method discussed in Section 2, These 
incremental values, along with the buffet onset information, were exercised 
in an extensive correlation study using known aircraft geometric properties 
and variations on the transonic simularity rules wherever these rules were 
deemed applicable. The additional information on current advanced or super- 
critical airfoil configurations Included in the correlation reflects the latest 
levels of technology. 

f 

Total aircraft drag is considered as being composed of two major parts, 
that which is Independent of lift . and drag due to lift or Induced 

drag (Cj)^). 


C 


D 



( 1 ) 


The minimum drag level occurs when drag due to lift is zero and is the. 
result of friction, compressibility, and miscellaneous items such as external 
stores. 



+ AC^ + AC 


( 2 ) 


Drag-due- to-lift is assumed to be composed of the incremental value of 
wing pressure drag and the theoretical level of Induced drag. 


- AC_ + C- VaR (3) 

°L ^ 

The terms and AC^p represent a departure from most classical drag 

buildup techniques and will now be defined in order to forgo any confusion 
regarding the above equation. The compressibility drag increment due to 
volume (ACdj,) is assumed to be composed of form drag which is a function of 
shape and volume effects upon viscous pressure levels due to Increases in 
local Mach numbers, and compressibility effects as local flows become sonic 
and shock waves form. In the supersonic speed range, form compressibility, 
and zero lift wave drag are assumed synonymous. 
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The wing pressure terra ) is the induced drag increment above that 

produced by the theoretical level resulting from It primarily is 

a function of the wing and Includes a combination of many effects (separation, 
lift dependent compressibility drag, span-wise flow, body effects, etc.) none 
of which are easily computed with the limited aircraft data available. The 
division of Induced drag into two parts, theoretical and an Incremental wing 
pressure term, provides the user an easy way of treating induced drag since 
the above lift dependent drag items are not usually known in a preliminary 
design environment. 

This section presents the final correlated values resulting from this 
emperical drag correlation technique, called the "Delta Method". These para- 
meters also make up the aircraft drag prediction computer program (EDET) 
detailed in reference 2* Included in this section is a discussion as to the 
use of these curves and their range of applicability. The actual data cor- 
relations are presented in Section 3. 

1.1 Design Lift Coefficient and Mach Number 

1 

The design lift coefficient (Ct ) of a configuration can be obtained 

^DES 

from the curves of figure 1 as a function of wing sweep, aspect ratio, effec- 
tive thickness and camber. The presentation of versus AR 

was derived from configurations with wing effective thickness ratio less 
than 0.065. The curve at the top of the page is primarily for subsonic/tran- 
sonic configurations where (t/c) is greater than 0.065 including current 
supercritical and advanced technology wings. The data correlation which pro- 
duced these curves is presented in Section 3, figures 30 and 31. The lift 
coefficient upon which the majority of the following drag computations are 
based is defined as the Increment from this design lift coefficient, given by 
the relationship 


AC - C - C, (4) 

^ ^ ^DESIGN 

After the design lift coefficient is determined by the above method, the 
two-dimensional drag divergence Mach numbers is next obtained from figures 2 
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and 3 at the design lift coefficient. Figure 2 is to be used for supersonic 
configurations and figure 3 dif ferenciates between conventional and the effects 
of the current supercritical and advanced wing technology for subsonic con- 
figurations. The data correlated to produce these relationships is presented 
in figures 32 and 33 of section 3. 

The two-dimensional drag divergence Mach number, is corrected to a three- 
dimensional value by making the corrections for sweep and aspect ratio effects 
as presented in figure 4 to derive the configuration's design Mach number. 

The values are the same as given by reference 33 and are applicable to super- 
sonic, conventional, and advanced sections. 

The design Mach number then determined by the relationship 

and the incremental Mach number upon which the majority of the following drag 
computation is based is determined by 

AM - - M^gg (6) 

1.2 Turbulent Skin Friction (Cp^) 

After the determination of the configuration's design and Mach number 
is completed, the next step is the computation of those component drag levels 
which are Independent of lift. The first of these is the determination of 
skin friction. 

Using a three-view or sketch of the configuration in question and esti- 
mated area progression curve, the wetted areas and reference lengths of the 
fuselage, wing, and tail surfaces are determined. The reference length for 
the fuselage is its overall length. This dimension for the wing and tail sur- 
faces is their individual exposed mean aerodynamic chords. The inlet capture 
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TM Should be rsflnsd in ths subject sre. progression to sssure . correct 
Physlcsl presentation for the clculetlon of oetted .res end fuselage fineness 
ratio. The reference altitude for drag estin«tlon or Reynolds nunber per foot 
together with the desired Mach nusd,er la selected by the user In computing C . 
This option la also offered In the computer code (EOET) In order that the pr^ 

gram may be used either to compute aircraft drag or to verify (or check) wind 
tunned results, 

lenethTl'°r*““‘°” reference 

lengths is then computed and the total configuration akin friction drag estl- 

mated by the expression 


^0 S 


® "■ (t) 


(7) 


where is the total number of configuration component parts. 

Cp. flatplate friction coefficient, Is determined from figure 5 which Is 


Cp is presented as a function of component 


extracted from references 34 and 38 

Reynolds number and transition location (,/f). Th. aero transition location 
la assumed for aircraft drag eatlmatlon. but a value of (x/f) muat be spec- 
ified with work concerning wind tunnel models where grit Is generally applied 
to assure a specified transition location. 

r.F. , Form factor. Is taken from figures 6 and 7 which are extracted 
from reference 34. Body form factor la presented In figure 6 as a function 
of fineness ratio (f/d). Care must be teken for the case of Internally 
mounted engines that Inlet capture area be Included to obtain the correct 
value of f/d. wing and tall form factor, are given by figure 7 as a function 
of effective thickness ratio. A distinction Is also made between conventional 
and advanced airfoil section. The reader 1 . referred to Appendix A for the 
erlvatlon of the Information presented In these two figures. 

Since the value of c^ obtained from figure 5 Is based on zero Mach num- 
ber. It ™.t be corrected for Mach effects. This Is accomplished by the In- 
c uslon of the term (Cp/Cp^^^) which Is the ratio of compressible to 
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incompressible friction drag obtained from figure 8 at the desired Mach number. 
These data are extracted from information as contained in reference 35. 

For most fighter configurations there is an additional level of miscel- 
laneous drag associated with items such as gun ports, gutters, antennas, blis- 
ters, and conopy shape. These items are Included in the flight test polars 
for each configuration presented in references 3 through 22. A drag buildup 
of these appendages could be made if a complete and detailed description of 
them were available, however this is not usually the case during preliminary 
design. A realistic overall drag level can successfully be estimated by a 
percentage increase in the level of computed friction drag. Figure 9 pre- 
sents a correction factor by which the miscellaneous drag estimation may be 
made. It is a result of the correlation of the computed friction drag of a 
representative sampling of the subject aircraft versus their recorded sub- 
sonic minimum drag levels. This correlation is presented in figure 34, Sec- 
tion 3 and the derivation of the correction factor is presented below. 

The value of average skin friction coefficient was computed from the 
relationship 


AVG . Cp ( comp uted ) 



( 8 ) 


where Cr. was computed from equation (7) , and is the total configure- 

Dp TOT 

tlon wetted area. The value of an actual average Cp for the configuration 

was then obtained directly from figure 9 and this value used to predict the 

actual subsonic friction drag level by the relationship 

( AVG. Cp actual \ 

AVG. computed") 

F / (9) 

“ 1.284 C (computed) 


where 1.284 is the slope of the curve as presented in figures 9 and 34. For 
the case of wind tunnel models, or when the level of configuration friction 
drag is already known, the above correction procedure is not necessary. 
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1.3 Compressibility Drag ) 

c 

Configuration zero lift compressibility drag (Wave) is assumed to be 
composed of three components - that due to wing, the fuselage, and wing/body 
interference effects* 


AC - AC WING + AC FUS AC INT 


( 10 ) 


Correlations for component compressibility drag are presented in the 
following discussion* 

1*3*1 Wing compressibility drag WING) * - Subsonic wing compressibility 

drag is presented in figure 10 as a function of effective thickness, camber, 
and AM. These curves are used for both conventional and supercritical and 
advanced technology airfoils. The data correlations to produce these relation- 
ships is given by figure 35. 

Figure 11 presents supersonic wing compressibility drag as a function of 
thickness, camber, aspect ratio, sweep, andAM* These curves are restricted 
to wings of effective thickness ratio less than 0.065 and aspect ratios less 
than five (5) as these were the outer bounds of the data from which the cor- 
relation was derived. The reader should note the change in correlation par- 
ameter between subsonic and supersonic configurations. The data correlation 
for figure 11 is given by figure 36 of Section 2. 

1*3.2 Fuselage compressibility drag (AC n^ FUS) . - To estimate properly the 
compressibility drag of a fuselage, the user of this technique begins with a 
preliminary estimation of the area distribution* Here, in contrast to the 
friction calculation, it is extremely important that this distribution not 
include the inlet capture area* The cross-sectional areas of horizontal and 
vertical tail surfaces should be included, however, in order to remain con- 
sistent with the methods used in developing the following data* The contri- 
bution of the empennage to the overall compressibility drag is usually small 
when compared to the whole but still should not be overlooked. 
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From the area distribution the user estimates the total length (£), 
maximum cross-sectional area (S-p-), and base area These parameters are 

configuration oriented and large variations can be expected from aircraft to 
aircraft. A consistent method of determining their values is, therefore, 
necessary to make this drag estimation procedure meaningful. Figure 12 pre- 
sents the method used in this study for fuselage geometry estimation. This 
method should be followed as closely as possible to maintain consistency 
throughout the drag computation. Area progression curves for the majority of 
the study aircraft are presented in reference 33* 

Subsonic fuselage compressibility drag is presented in figure 13 as a 
function of the geometry parameters and Mach number. It should be noted that 
fuselage and interference drag computations are based on the free stream Mach 
number, while all other drag items are based on the design Mach and are 
correlated against AM. The value of obtained from figure 13 is based on 
maximum fuselage area (5,^) and is converted to the reference wing area by 
the relationship. 



Figure 14 gives supersonic fuselage compressibility drag as a function 
of body fineness ratio, base to maximum area ratio, and Mach number. Equa- 
tion (11) must be used to convert the value of obtained from this figure 
to the correct area relationship also. The data correlations which produced 
figures 13 and 14 are presented in Section 3, figures 37 and 38, 

1.3.3 Interference drag (ACp ^ INT) . - Wing/fuselage Interference drag, as 
derived for this method, is presented in figure 15 as a function of body 
diameter to wing span ratio, taper ratio, and wing sweep. The assumption 
la made that interference drag is zero below Mach number equal to 1.0. The 
correlation of data which produced this relationship is shown on figure 39. 

1 . 3.4 Total compressibility drag ) « ” Th® values of AC^^WING, ACd^FUS, 

and ACj, INT are now combined to produce the total configuration compressibility 
drag by Sse of Equation (5). The resulting values, when added to the previously 
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computed level of frlctlou dreg, will produce . conflguretion dreg level which 
Is Independent of lift throughout the desired Much renge. 


1«4 Miscellaneous Drag 

Allowance is made within the coding of this procedure for the addition 
of a miscellaneous drag level (AC^,) so that the effects of a configuration 
change or an external store may be Included, Such drag levels are computed 
external to this procedure and input as a function of Mach number. The 

present coding of EDET allows for a total of 10 separate miscellaneous drag 
items. 


1*5 Induced Drag 

Lift dependent drag is predicted as the sum of the theoretical Induced 
drag (Cj^ /ttAR) and a derived term called wing pressure drag (ACpp) . The 
theoretical value is used to define the primary variation of induced drag 
with lift coefficient and assumes an Oswald's efficiency factor (e) equal 
to unity. The wing pressure term is, therefore, composed of all lift de- 
pendent drag items over and above that of theory. This division provides 
the user with an easy method for treating induced drag in a preliminary 
design atmosphere where not much is usually known concerning a configuration's 
drag characteristics with variations in angle of attack. 

1.6 Wing Pressure Drag ^Cjjp) 

The remaining lift dependent drag, is Included by this technique in the 
wing pressure drag term, ACop. Figures 16 through 20 present the subsonic 
pressure drag as a function of aspect ratio, thickness, camber, AM, and AC . 
The curves are entered at the desired value of M and interpolated between^t 
the required value of AR for the range of AC^^ under consideration. 

These data apply to current supercritical and advanced technology wing design 
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as well as that of conventional subsonic sections. The correlation of data 
for these curves is presented in Section 3, figures 40 through 48. 

Figures 21 through 27 presents a similar correlation of wing pressure 
drag for supersonic airfoils as a function of the same variables as the sub*- 
sonic data. Interpolation between curves is handled in the same manner. 

These curves, however, are restricted to an effective wing thickness of 0.065 
and values of AR less than 2.0 which are the outer bounds of the 

basic data from which they were generated. Data correlation for these curves 
is presented in figures 49 through 57. 

1.7 Total Aircraft Drag (C^) 

The drag items computed under Sections 1.2 through 1.6 are now combined 
with the theoretical induced drag levels by the relationship 

C - C„ + AC_ + AC- + AC + C ^/ttAR (12) 

D Dj. D^ D °p ^ 

The resulting values, then, represent the estimated total configuration 
flaps up drag at the Mach number and lift coefficient selected by the user. 


1.8 Buffet Onset (^Lg^O,^ 

The drag estimation technique described in the preceeding sections may 
also be used for predicting the lift coefficient for buffet onset 
To this end the buffet onset Information for the subject aircraft were 
evaluated in much the same manner as was the component drag levels. Figure 28 
is used for predicting Cl^ ^ as a function of thickness, aspect ratio, camber, 

sweep, ^ES’ 
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It can be seen that buffet onset is relatively independent of wing 
thickness at low speeds but, as the design Mach is approached, it begins to 
show an ever increasing variation. The end result of this correlation 
indicates that, once the design Mach number is passed, the thicker wings 
tend to enter into buffet at values of lift lower than while the 

thinner wings common to present day fighter configurations exhibit an 
increasing margin in load factor from the design lift within the Mach number 
range of the basic data. 

It should be pointed out that this correlation does not include the 
effects of maneuvering devices such as flaps or slats but is representative * 
of the clean wing only. The data correlation that resulted in this relation- 
ship is presented in figure 58, Section 3. 
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Figure 5. - Variation of flat plate incompressible turbulent skin friction coefficient with 
Reynolds number. 











Figure 7. - Wing section form factors 






Figure 8. - Compressibility correction to skin friction coefficient 




























Figure 15. - Wing/body zero lift interference drag 
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Figure 19. - Subsonic wing pressure drag, AR ' ■ 4.0. 
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Figure 20. - Concluded. 
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2. DATA REDUCTION TECHNIQUE 


\ 


The aircraft drag information, as presented in references 3 through 22, 
was analyzed in accordance with procedures as suggested by the Lockheed 
developed Delta ^^ethod of references 1 and 34, This method suggests that 
when incremental values of component drag are compared in relation to the 
increment (A) in Mach number and lift coefficient from the configuration's 
aerodynamic optimum values (Mpgg and ^ collapsing of the data is 

produced when transonic similarity parameters are utilized and a reasonable 
correlation of the data is therefore possible. The method of drag analysis 
is presented graphically in figure 29 and outlined in the following discussion'. 

2.1 Design Lift Coefficient and Mach Numbers 

For each aircraft under consideration, the reference drag data were 
used to construct curves of Mach number times lift-to-drag ratio, M(L/D) , 
for a given value of lift coefficient from the relationship 

M(L/D) - M(Cj^/Cj^) (13) 

These quantities were then plotted versus Mach number and the maximum 

value of M(L/D) determined for each value of C , Ninety-nine percent of 

L ^ 

this maximum value of was then determined and the Mach number at which 

it occurred noted for each C . These quantities were then plotted as 

T L 

0.99Mhr) versus C, and M. The peak of the curve of 0.99Mf;7) versus 
D max L D max 

was assumed to define the design lift coefficient and the corre- 
sponding Mach number, the design Mach should be noted that these 

conditions define the aerodynamic optimums only and’ may or may not reflect an 
optimum flight cruise condition when engine SFC is introduced into the maxi- 
mum range factor M(-^)/SFC. Examples of this procedure are given by refer- 
ence 33 for each of the study aircraft. 
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2,2 Drag Divergence Mach Number 


The three-dimensional Mach number for drag divergence is defined as the 
point where the rate of change of drag with Mach number, (dCL /dM) , is equal to 
0.10. These values of Mach number were noted from curves of total drag coef- 
ficient versus Mach number for constant values of lift. Corrections to a two- 
dimensional value were made using the method contained in reference 34 which 
is a function of wing aspect ratio and quarter chord sweep angle. Examples of 

the resulting parameter can be seen for each aircraft in Section 3 of this 
report. 


2 . 3 Drag Breakdown 

The drag of any aircraft configuration can be broken down into two basic 
components - that which is independent of lift, and that which is the direct 
result of lift generation or angle of attack. The first of these (lift inde- 
pendent drag) can be thought of basically as the aircraft's resistance to 
movement and constitutes the minimum drag level of the configuration. The 
contributors to minimum drag go by many names but generally consist of 
friction drag, form or pressure drag, and compressibility effects. Friction 
drag is a function of sire. Reynolds Number, and surface roughness and is 
associated with shear stresses in the configuration boundary layer. Form or 
pressure drag is affected by shape or volume effects such as thickness, fine- 
ness ratio, and contour slopes and is associated with increases in viscous 
pressure levels due to increases in local Mach numbers. Commpresslbillty 
effects become apparent with Increased speed as local flows become sonic and 
the formation of shock waves occur. In the supersonic regime, form, com- 
pressibility, and wave drag are synonymous. For the purposes of this study, 
form and compressibility drag have been lumped together into one term called 
compressibility drag due to volume (AC^^). The minimum drag level, therefore, 
is assumed as composed of only two parts - friction drag which is easily cal- 
culated. and compressibility drag. The procedure by which these two drag 
levels are determined will be discussed in succeeding paragraphs. 
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The lift dependent (induced) drag level is a function primarily of wing 
size and shape and is associated with the spanwise load distribution generated 
and its variation with angle of attack and speed. Depending on the configura- 
tion, there is also a certain amount of compressibility or wave drag due to 

produced as angle of attack is increased. This compressibility drag 
increment is usually small, however, when compared to the total induced drag 
level but is included as a contributor to lift dependent drag. 

For the purpose of this study, induced drag is divided into two compon- 
ents - the theoretical value (C /ttAR) which assumes an efficiency factor (e) 

la 

of 1.0 and an elliptical spanwise load distribution, and a wing pressure term 
(ACjjp) . Compressibility effects due to angle of attack are Included in the 
ACpp term. The procedure by which these two drag values are determined is 
presented in the succeeding discussion. 

For each aircraft under consideration, tabulations of total configura- 
tion drag were made as a function of Mach number and lift coefficient. These 
tabulations did not include the effects of wing flaps, landing gear, or 
external stores and were, therefore, representative of a clean configuration. 

As the first step in analysis, the theoretical induced drag term was 
removed by the relationship 

c; - C,. - C-^/ttAR (14) 

total 

where the primed value represents an Interim calculation. 

A low speed Mach number, usually M - 0.60, approximately AM - -0.30 
below the design Mach was then chosen and the resulting polan shape obtained 
from Equation (lA) plotted. The minimum value of drag obtained from this 
operation was designated Cq— or minimum parasite drag. It was assumed that 
this value was composed entirely of friction drag at the Mach number chosen; 
although it is realized that in actuality, this value usually includes some 
pressure drag. Skin friction drag at other Mach numbers was then computed 
using the relationship 
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^ ^INC “ 


P min I (C„/C^ ' ) 

F ^INC M » 0.60 


( 15 ) 


Values of are presented in figure 8 and are based on information 

contained in references 34 and 35. 

These computed values of skin friction drag were then removed from the 
results of Equation (14.) 



(16) 


The resulting lift/drag relationships were plotted for each Mach number and the 
minimum drag quantity obtained was then identified as the compressibility or 

zero-lift wave drag (AC^^^) of the configuration inclusive of wing. body, and 
wing /body interference. 

Wing pressure drag was then computed from the relation 


and its variation with Mach number and lift coefficient noted. As mentioned 
earlier, the term AC^p Includes the combined effects of flow separation, com- 
pressibility drag due to lift, and trim drag. Any attempt to further break the 
term down into its components was judged to be beyond the scopd of this 

study and would require levels of detail outside the Intended preliminary 

design application. Values of AC are presented for each subject aircraft 
in reference 33. 
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2.4 Component Compressibility Drag 


The values of zero-lift wave drag obtained in the above discussion are 
those of the entire configuration. The three main contributors to total con- 
figuration wave drag are the wing, fuselage, and wing/f uselage interference. 

AC - (AC )WING + (AC )FUS. + )INT. (18) 

^C “c C 

To separate out these component parts, the following procedure was used; 

• The fuselage area distribution, minus capture area, was determined 
from information as contained in references 3 through 21. Horizontal 
and vertical tail cross-sectional areas were Included as basic fuse- 
lage area since their contribution to the total wave drag is usually 
small but should not be neglected. 

• This area distribution was then input to the NASA Wave Drag Program, 
reference 36, as an equivalent body of circular cross section. The 
wing was then added to this area distribution and the zero-lift wave 
drag of the resulting configuration was computed at supersonic speeds 

• From the output of the Wave Drag Program, incremental values of 
compressibility drag due to the fuselage, wing and wing/f uselage 
interference were obtained and the percentage of total wave drag 
for these items de termined versus Mach number. 

0 Xhese percentages were then applied to the values of ACj)^ obtained 
from the results of Equation (16) to produce incremental values of 
supersonic compressibility drag due to fuselage, wing, and fuselage/ 
wing Interference. 

Subsonic and transonic fuselage compressibility drag and drag divergent 
Mach number were estimated using results obtained from wind tunnel tests of 
bodies of various fineness ratio, reference 34. These values were then sub- 
tracted from the total value of ACq derived from the basic data within this 

c 

Mach range. The remaining quantity was assumed to be wing compressibility 
drag. Values obtained from this procedure were then faired into existing 
supersonic drag levels as computed in the preceeding paragraphs. For most 
of the aircraft considered, the fuselage drag divergence occurred at a Mach 
number sufficiently high enough to allow the values of computed to be 

due entirely to the wing at subsonic speeds. Interference drag was assumed 
to be zero below Mach 1.0. 
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3. DATA CuRR£LATI0i4 


Sections 1 and 2 have i>resented the Delta Method for estimating total 
aircraft drag^ an overview of final results of the correlation effort, and 
the method to which the test data were broken down. This section presents 
the individual aircraft /wind tunnel model data points for each correlated 
parameter along with the rationale for selection of correlation variables. 

The original development of the Delta Method Drag correlation technique, 
described in reference 1, was directed toward determination of drag correla- 
tion techniques which would permit preliminary design evaluation of the 
potential of transonic aircraft incorporating advanced or supercritical air- 
foils, This current work expands the data base to include conventional and 
advanced airfoils operating in the transonic flight regime and conventional 
airfoils operating in the supersonic flight regime. At the present time, the 
NASA Air Force TACT aircraft represents the only available /known data where 
advanced airfoil design has been incorporated into a supersonic aircraft. 

These flight data became available too late to be incorporated into this report- 
ing; however, incorporation at a later date would represent a valuable addition 
to the Delta Method. 


In developing the various component drag correlations, extensive guidance 
was derived from previous work by McDevltt (reference 37), DATCOM 
(reference 35), and RAS Data Sheets, The effect of external stores, tip 
tanks, flaps, cruise droop, and other camber varing devices were not considered. 
A clean wing was assumed at all times. 


The basic geometric properties of all 19 study aircraft are presented 
in table 1 along with 15 current supercritical or advanced technology experi- 
mental wing designs. Included are the values of and for 

the study aircraft calculated by the methods discussed in Section 2 and shown 
in reference 33, The plotted symbol by which each configuration is recog- 
nized is also presented. The flagged symbols indicate wind tunnel model data 
extracted from reference 1. On the correlation figures to follow, the darkened 
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symbols represent the advanced or supercritical model designs on reference !• 
Correlation parameters used in the analysis are tabulated in table 2. 

It has been assumed that the drag polars presented for all aircraft or 
models are trimmed, and at the low lift coefficients of interest, trim drag 
is negligible except at supersonic speeds where it is included in the wing 
pressure drag term. Considering the diversity of designs, i.e., supersonic 
fighters, subsonic attack aircraft, surveillance aircraft and transport 
models varying contractor sources of data, differences in bookkeeping methods, 
power effects, static margin, and varying test Reynolds numbers, the component 
data correlation is considered to be quite adequate for the preliminary design 
of both subsonic and supersonic aircraft, 

3.1 C and 
DES 

Design lift coefficient and design Mach number for each configuration 

were defined as the subsonic C and M at which 0,99 M (L/D) occurs. 

^ MAX 

(Note that these are aerodynamic optimums and may or may not reflect an 
optimum flight cruise condition where engine SFC must be Introduced in the 
maximum range factor M (L/D)/SFC). 

3.1.1 design. - Subsonic wing design techniques consider the variables 
of aspect ratio, taper ratio, sweep, thickness ratio, camber, and twist in 
"designing to" optimum pressure distributions in deriving C and M design. 

It was premised that for these subsonic aircraft, an elliptical span load at 
the design conditions would be attained, therefore, twist and taper need not 
be considered as unique correlating parameters. Further, subsonic wing design 
practices rely heavily on a high degree of two dimensional flow for accurate 
prediction of design, hence some minimum aspect ratio could be anticipated. 
The correlation on figure 30 reflects the relationship of the remaining primary 
design variable to design for both advanced and conventional wings. The 
aircraft represented in this correlation have been tailored for best cruise 
performance and compromised only by subsonic cruise speed and take off and 
landing considerations. 
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With Capture Area roM>vad for Internally Bounced wiglnea 

Based on 235 In. diem fuselage excluding vheel fairings and center engine inatallacion 




Optimum cruise efficiency for supersonic tectical aircraft geometries 
may be compromised by maneuver, takeoff and landing, ceiling, external store 
carriage, and maximum speed. The high dash speed requirement drives those 
designs to low aspect ratio-highly swept-thin wings with resulting aerodynamic 
best cruise efficiency occurring at high subsonic speeds, sometimes higher 
than the best cruise speed for the engine. For this class of aircraft, 
design (see figure 31), correlates with aspect ratio times thickness ratio 
to the one-third power. The change over in design relating parameters 
occurs at approximately an aspect ratio of five. (Note the T-2B and A-6 on 
figures 30 and 31). 

3.1.2 M design. - The original Delta Method derivation (reference 1 and 34) 
shows that M design and M divergence, i.e., d C^/dM = 0.10, are synonymous 
at the design lift coefficient. M design for the subsonic configurations 
(see Figure 32), when adjusted for the three-dimensional effects of aspect 
ratio and sweep, follow the correlating parameters suggested in the Transonic 
Similarity Rules (reference 37). The increased drag divergence Mach number 
for those configurations employing advanced sections is evidenced. For the 
supersonic aircraft designs, i.e., low aspect ratio-highly swept-thin wings, 
Mach design corrected for three-dimensional effects, correlates more favorably 
as a direct function of wing thickness ratio. (See figure 33.) 

Referring back to figure 32, this correlation Implies that as the wing 
thickness is increased above a thickness ratio of approximately 20 percent, 
little benefit in drag divergence Mach number is evidenced for the advanced 
airfoils. This logic is borne out in the limit case of a cylinder where 

- 1.0. There would be no way of dictating differences in upper and 
c 

lower surface pressure distribution and no distributions in section charac- 
teristics. Until advanced sections of greater thickness ratio are developed, 
this diminishing return of advanced over conventional sections appears valid 
for preliminary design purposes. As the thickness ratio is reduced, the 
spread in divergence Mach number between the two technologies is increased. 
Again, in the limit case with thickness ratio approaching zero, no mechanism 
for prescribing desired upper and lower surface pressure distribution exists. 
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Hence relations of advanced airfoil section drag divergence Mach number 
should not be extended below the thickness for which data are shown- This 
airfoil technology level difference for lower thickness ratios does lend 
emphasis to consideration for advanced section technology application to 
supersonic fighters for a potential improvement in cruise performance. 

3-2 Friction Drag Correction Factor 

The Delta Method includes a procedure for determination of the subsonic 
friction drag on bodies, wings, tails, and nacelles- This theoretic friction 
estimate must be increased by approximately six percent to produce friction 
drag levels that agree with flight test results of transport and bomber type 
aircraft. For attack and fighter type aircraft which might have bubble 
canopies, weapons, pylons, antenna in flight refueling probes, etc,, this 
factor is much larger, A correlation of average skin friction computed by the 
theoretic procedure is compared with the actual average skin friction 
coefficient in figure 34 to derive a simple and more realistic correction 
factor for fighter/attack type aircraft friction drag. This factor is 
28.4 percent. When using the Delta Method, good engineering judgment should 
be exercised in choosing the procedures for friction drag determination most 
appropriate to the concept under consideration. The correlation for the 
friction drag correction factor presented in figure 34 is the same as given 
by figure 9. 


3.3 Compressibility Drag 

Compressibility, wave or pressure drag, i.e., all synonymous terms, 
generally is separated as drag due to volume, i.e., zero lift, and drag due 
to lift Compressibility drag due to volume for this correlation has been 
broken down into the contribution of the wing, the fuselage, and a wing body 
interference contribution (see figure 29, Section 2.0). The compressibility 
drag due to lift for both the wing and fuselage are Included in the pressure 
drag correlation explained in section 3,3,4. 
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j.3.1 Compressibility drag due to wing volume. - The only transonic and 
supersonic linear wave drag theories available generally prescribe a syste- 
matic series of uncambered wings having specific section characteristics, 
i.e., diamond, double wedge, biconvex, etc. Since the wave drag required 
for this method was for real airfoils and no systematic sets of wing 
variables existed, selection of the correlating parameters suggested by the 
transonic similarity laws and supersonic drag theory were used on a best fit 
of data approach. Transonic similarity laws suggest normalizing WING 

by (t/c)^^^ which was done. Based on previous correlations, reference 1, 
the wing compressibility or wave drag was also normalized by a camber factor 
of (1 + — /lO). In the subsonic/transonic speed regime, compressibility or 
wave drag for the wings of advanced or conventional sections were found to 
correlate best as a function of (t/c)^''^ (see figure 35). In the supersonic 
speed regime, AR tan A proved to be the correlation parameter that pro- 
duced the best data fit for all airplanes (see figure 36). The F-104 and 
F-106, having a major wing geometry difference in sweep and taper only, were 
the drivers in selecting this correlating parameter. The final data fairings 
of figures 35 and 36 are summarized in figures 10 and 11. 

3.3.2 Compressibility drag due to fuselage volume . - Since fuselage drag 
characteristics are not dependent on wing design characteristics, fuselage 
compressibility or wave drag has been correlated at discrete Mach numbers 
rather than at a AM from wing design Mach number. The procedures followed 
in separating wing and fuselage compressibility or wave drag due to volume 
from the total aircraft minimum drag have been identified graphically on 
figure 29, Body or fuselage wave drag normally is divided into three compo- 
nents, i.e., forebody, aftbody, and midbody interference drag. Linear 
theory provides a guide to data correlation through the relationship 
C (2/d)^ - (f ) M^ - H aid ) . Initial correlation attempts to separate 
total fuselage drag into these components proved fruitless, partially due 
to the extreme difficulty in defining the extent of forebody, midbody, and 
aftbody from the area progression curve. Observing the many area progression 
curves, those parameters that appeared to be consistently definable for all 
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configurations were the overall length, the maximuni diameter and a base area. 

The base area was not interpreted - soley as a base drag contributor - it 
Is also an Indicator afterbody closure, i.e., the higher the base area the 
less severe an aftbody closure and the lower the aftbody wave drag contri- 
bution. From these observations and assumptions, fuselage wave drag 
(ACj) FUS.) was multiplied by the square of the length to diameter ratio, 
and correlated as a function of one plus the ratio of base area to maximum 
frontal area (see figures 37 and 38). These data fairings are used in 
figures 13 and 14 of the methodology. 

333 Wing body interference dra g. - Wing-body interference drag was 

determined from the Lockheed modified NASA v;ave drag program of reference 36. 

Interference drag is defined as the difference in total aircraft equivalent 

body wave drag and the sum of the theoretic component, wing and fuselage 

drags taken separately. The wing-body Interference drag curve on figure 39 

exhibits both positive and negative interference. Wing geometries are the 

major quantifiable contributor to differences in this interference parameter. 

Various combinations of wing geometries were investigated to achieve the 

best fit of data. Interference wave drag is normalized by (1 - X) times the 

Cos A and correlated as a function of wing span to body diameter. The 
c/4 

solid curves of figure 39 show the final data fairings which are summarized 
in figure 15. 

3.3.4 Wing pressure drag due to lift. — The Delta Method defines the pressure 
drag term as additional lift related drag over and above the classical Induced 
drag represented by C ^/irAR. This pressure drag results from flow separation 
and shock- induced losses, and in the supersonic speed regime, the drag term 
ACpp is also comprised of some fuselage lift. 

Both transonic and supersonic linear theory predict drag due to lift, 

fc ^/C Vo be a function of leading edge suction which suggested that wing 
\ b/ 

geometries that effect leading edge suction would produce the best correlation. 

Transonic similarity laws gave the clue as to the correlating parameters to 

1 /3 

use. Best data fit was obtained by first normalizing AC^p by (t/c) times 

the camber term (1 + ~ /lO), and correlating versus AR (t/c) . Correlation 

c 
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of this term is presented on figures 40 through 48 for subsonic con- 

figurations and on figures 49 through 57 For supersonic configurations. 
The solid curves showing the final data fairings are used in figures 16 
through 27 of the Delta Method. Pressure drag due to lift correlation is 
continuous throughout the Mach number range and no difference exists be- 
tween conventional and advanced airfoil pressure drags when compared for 
identical geometries. 


3.4 Buffet Onset 

The buffet onset data were available for only a limited number of 
configurations and are restricted to wings of conventional sections with 
the exception of wind tunnel buffet measurement on the T-2C supercritical 
wind tunnel model. Buffet onset is very difficult to define due to diffi- 
culty in determining intensity and the lack of common threshold in- 
tensity definition. Considering this, the correlations shown on 
figure 58 are considered excellent. Until more advanced section buffet data 
becomes available, it may be premised that this correlation is appropriate 
for both conventional and advanced or supercritical wings. Data are shown 
as a function of wing sweep, thickness, aspect ratio, camber, ^L])es» 

AM. The Mach number range for this analysis is -0.6 < AM < +0.20. The 
solid curves represent the final data fairings which are summarized in 
figure 28. 
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Figure 32. - Mp2-D correlation, subsonic configurations. 
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Figure 32. - Concluded. 
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Figure 33. - Concluded. 
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Figure 36. - Continued. 
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Figure 39. - Continued 
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Figure 41. - Concluded. 
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Figure 42. - Subsonic ACj^p correlation ACl “ 
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Figure 43. - Concluded. 
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Figure 44. - Concluded. 


99 




100 




AR (Vc) 


1/3 


Figure 45. - Concluded. 
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i-igure 47. - Concluded, 
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Figure 48. - Concluded. 
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Figure 50 
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Figure 50, - Concluded* 
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Figure 51. - Continued. 
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Figure 52, - Supersonic wing pressure drag correlation, AC, « -0,05, 
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Figure 52. - Continued 
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Figure 54. - Supersonic wing pressure drag correlation, - +0.05. 
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Figure 54, - Continued 
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Figure 55. - Continued. 
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Figure 55. - Concluded. 
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Figure 56. - Concluded. 
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Figure 57, - Supersonic wing pressure drag correlation, = +0.30. 
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4. EXAMPLE AIRCRAFT DRAG BUILDUP 


The Delta Method drag prediction technique, as presented In Section 1, 
is used in this section to compute the drag polars of several of the subject 
aircraft at selected Mach numbers. These examples are then compared to the 
basic data for the subject aircraft as presented in references 3 through 21 
and discrepancies noted. Included also is a comparison of certain drag 
levels for an advanced supersonic configuration which was not Included in the 
original stable of subject aircraft. This method has also been converted 
into a computer code named Emperlcal Drag Estimates Techniques (EDET) . This 
code is compatible with the NASA-Ames computer facilities and is detailed in 
reference 2, 


4,1 Design Lift Coefficient and Mach Number 

The first step in the Delta Method is the determination of the design 
lift coefficient and Mach number of the configuration in question. This is 
determined from the values of aspect ratio, thickness, sweep, and camber of 
the configuration being analyzed. 

The calculation for is given in table 3. The curves of figure 1 

UtiD 1 / 

are entered at the appropriate value of AR (t/c) for supersonic and AR for 
subsonic designs and the values of noted. For the subsonic configu- 

rations, the value taken from the curve must be multiplied by the expression 
cos Ac/4 (1 + ^q'^Vn/aR. Also included in table 3 are the values of 

derived from the basic data and used in figures 30 and 31 to correlate 

It can be seen that the deviation between these values and those derived from 
figure 1 are (with three exceptions) within a band of ACj^ - ±0.05 which is 
considered acceptable. 
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TABLE 3, - DESIGN LIFT COEFFICIENT CALCULATION 
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The calculation for is given by table 4. The curves of figures 2^^ 

and 3 are entered at the appropriate value of ^/c for supersonic and ( /c) 
for subsonic designs and the corresponding value of Md 2 _o noted. For subsonic 
configurations the distinction must be made between advanced or conventional 
airfoil sections. The above value of two-dimensional drag divergent Mach is 
then corrected for sweep and aspect ratio using figure 4 and equation (5) . The 
resulting three-dimensional drag divergent Mach number is, then, the design 
Mach (Mjjgg) . It can be seen from table 4 that the deviation between these 
calculated values and those derived from the basic data are (with two excep- 
tions) within a band of M - ±0.05 which is considered acceptable. 


4.2 Friction Drag 


Once the configuration design lift and Mach number have been determined, 
the next step is to calculate those drag values which are independent of lift 
and which make up the minimum drag level. The first of these items is fric- 

tion drag. 

Friction drag is calculated as discussed in Section 1.2 using equa- 
tion (7). Two aircraft (AA-F and RA-5C) are chosen as examples to demon- 
strate the method. Component wetted areas and reference lengths are either 
obtained from the geometric data of these aircraft or measured from the 
three-views as presented in references 3 through 21. The computation is 
shown in the following paragraphs. First, the required geometric parameters 
are computed as follows: 



A-4F 

KA-5C 


Conponent 

^WET ft^ (m^) 

£ « ft ft (m) 

^WET/®REF 

SreT/®RET ft* (b*) 

£-ft (m) 

^WET/^REF 

Wing 

430.0 (39.95) 

9.48 ( 2.89) 

1.65 

1144.08 (106.29) 

13.95 ( 4.25) 

1.63 

HOTZ* Tall 

84.40 ( 7.84) 

4.66 ( 1.42) 

0.32 

388.72 ( 36.11) 

9.73 ( 2.97) 

0.56 

Vert. Tall 

49.95 ( 4,64) 

7.38 ( 2.25) 

0.19 

235.33 ( 21.86) 

8.34 ( 2.54) 

0.34 

Fuselage 

547.65 (50.88) 

40.70 (12.41) 

2.11 

1474.00 (136.91) 

73,30 (22.34) 

2.11 
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TABLE 4. DESIGN MACH NUMBER CALCULATION 
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Reference lengths for the wing and tail surfaces are chosen as the exposed 
aerodynamic chord lengths. The actual fuselage length Is used for the 
fuselage. 

Form factors are computed from figures 6 and 7 as follows: 


Component 

A - 4F 

RA - 5C 

Wing 

Horz. Tail 
Vert. Tail 
Fuselage 

*^/c « 0.0748 FF - 1.255 
*^/c - 0.055 FT = 1.185 
*^/c » 0.055 FT - 1.185 
£/d - 7.81 FF - 1.191 

*^/c - 0.050 FF - 1.165 

*^/c » 0.040 FF - 1.13 

*^/c - 0.040 FF - 1.13 

£/d - 8.52 FF - 1.155 


Fuselage fineness ratio (^/d) is computed with the-, effects of inlet capture 
area Included in order to obtain a true representation of fuselage geometry. 
Conventional airfoil sections are assumed for both aircraft. 


A reference altitude of 36,152 feet (11,019 meters) end a Mach number of 
0.60 Is selected for computation. A value of compressible to incompressible 
drag (C„/CrT„„) is taken from figure 8 as 0.966. A Reynolds Number per foot 

r XN* 


of 1.381 X 10^ is obtained from the atmospheric tables. For the case of wind 
tunnel model evaluation, a desired Reynolds Number and transition location in 
percent length must be specified. Here it is zero. 


Calculation of the friction drag is then obtained as follows: 


(Equation 7) I 

A-4F 

RA-5C 





®WET /S Vf/FFI 


r 

!5!l(f!_W(FF) 

Component 

RN X 10® 

(Fig. 5) 


RN X 10® 

S 

(Fig. 5) 

®REF'^FInc/ 

-^Sf 

Wing 

13.10 

0.00281 

0.00562 

19,27 

0.00264 

0.00484 

Horz. Tail 

6. 44 

0.00319 

0.00117 

13.44 

0.00280 

0.00171 

Vert. Tail 

10.19 

0.00294 

0,00064 

11,52 

0.00287 

0.00110 

Fuselage 

56.21 

0.00225 

0.00585 

95.29 

0.00219 

0.00516 

Total Cp 

F 



0.01328 



0.01281 
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These values are then corrected for configuration effects by the use of 
figure 9 as follows: (Equation (9)) 



A-4F 

RA-5C 

(Computed) 

0.01328 

0.01281 

(T) ■ From Figure 9 

1.284* 

1.284* 

” '^D (Computed) 

F F 

0.01705 

0.01645 

*Thls is the slope of the curve given by Figure 9. 



4.3 Compressibility Drag ^ " 

The second major contributor to zero-lift drag is that level due to 
compressibility effects. As discussed in Section 3, compressibility drag 
is estimated from the sum o7 the three major contributing parts. For this 
calculation the same two aircraft are chosen as were for the friction drag 
computation. Compressibility drag is calculated at Mach numbers of 0.9 and 
2.0 for the RA-5C and at 0.8 Mach only for the A-4F. 


^•3*1 Fuselage compressibility drag . - Compute the required fuselage param- 
eters as follows: 


Using Tables 1 and 2: 

A-4F 

RA-5C 


1.05 

1.27 


69.06 

93.32 

^TT^^REF " © 

0.073 

0,065 
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Fuselage compressibility Is then computed as shown. 


Using Figures 13 and 14 

> 

1 

RA - 

5C 


0.80 

0.90 

2.00 

C_ 4)^ (Fig. 13) 

- S4^/d)V(i) 

0 

0 

0.680 

0.00729 


- 

FUS - C X © 
C ^ 

0 

0.00047 


- 

(~y (Fig. 14) 

UjT a . 


_ 


18.25 

0.19556 

Ar rus - tL * © 
c ^ 


- 


0.01271 


It should be noted that the calculation for and fuselage fineness ratio 
(^/d) do not include the effects of inlet capture area. 


4.3.2 Wing compressibility drag. -- The following parameters are tabulated 
either from table 2 or 4: 



A - 4F 

RA - 5C 

Mjjgg From Table 4 

(1 + - © 

AR Tan A 

Lti 

0.875 

0.01430 

2.54 

0.178 

0.889 

0.00679 

3.49 

0.136 
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Wing compressibility drag Is then calculated as follows: 




RA 

5C 



- 0.90 

- 2.00 

CD 

AM - 

-0.075 

0.011 

1.111 

AC / 0 (Fig. 10) 

0.037 

0.220 


C 




AC WING - Fig. 10 X 0 

Q 

0.00053 

0.00149 

- 

AC / 0 (Fig. 11) 

-■ 

— 

0.675 

C 




AC WING - Fig. 11 X 0 

- 

- 

0.00458 

C 





4.3.3 Wing/body Interference drag . ~ Interference drag is assumed to occur 
only above Mach 1,0; therefore, only the supersonic calculation for the 
RA-5C is required. This is accomplished as follows: 


- 2.0 

00 

d/j^ - 0.1485 

(1 - \) CosA^^^ - 0.6427 - 0 

from Figure 15; (INT) x 0- 0.00048 

ACjj (INT) - 0.00048/^0 - 0.00075 

c 

9 

4.3.4 Total compressibility drag . - The three components of compressibility 
drag are then summed to obtain the total compressibility drag using equa- 
tion 10. These are compared with the original values of AC^^, obtained from 
the basic data (Equation 10). 









A - 4F 

RA 

- 5C 

M® - 0.80 

Mco “ 0.90 

M® - 2.00 

ACdj, fus 

0 

0.00047 

0.01271 

ACDj, wing 

0.00053 

0.00149 

0.00458 

ACdj, int. 

0 

0 

0.00075 

TOTAL 

0.00053 

0.00196 

0.01804 

ACd^ From Data 

0.00130 

0.00190 

0.01530 

Difference 

-0.00077 

+0.00006 

+0.00274 


Values of computed aircraft compressibility drag versus those obtained 
from the basic data of references 3 through 21 are compared in figures 59 and 
60 for some of the subject aircraft. It should be noted that the trends 
shown follow very well the reduced data, which is critical for methods used in 
design work. The compressibility drag of a more advanced technology fighter 
configuration, which is not included in the subject aircraft, was computed 
as 0.00194 at Mach 0.90 and 0.02493 at Mach 1.6 by this procedure. Analysis of 
available drag data on the aircraft reveal thse values to be 0.0017 and 0.2610 
respectively. 


4.4 Total Minimum Drag 

The total drag level for the configuration which is independent of lift 
(minimum) is now determined from a summation of the component drags computed 
in the preceedlng paragraphs. This calculation is as follows: (Equations (2) 

and (15)). 
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A-4F 

ra« 

-5C 


M®- 0.80 

Moo - 0.90 


Cp 0 M^- 0.60 - 

V^^INC @ o.eo-Q 

0.01705 

0.966 

0.01637 

0.966 

0.01637 

0.966 


0.949 

0.939 

0.768 

■ @ / ® - @f’a5)i 

0.982 

0.972 

0.795 

-“t, - ci » ® - ® 

F F 

0.01675 

0.01591 

0.01301 

^Cdc - © 

0.00053 

0.00196 

0.01804 

ACp - © 

0.00100* 

0 

0 

^Dmin -© + © + © 

0.01828 

0.01787 

0.03105 

basic data 

0.02040 

0.01599 

0.02800 

i 

Difference 

-0.00212 

0.00188 

0.00305 

Error 

-10% 

11.7% 

11% 


*The miscellaneous drag item (^j^) for the A-4F in the above calculation 
represents an estimation of the drag for a centerline rack and IFR probe 
which are included in the basic drag data as presented in reference 6. 


Agreement is acceptable for the two aircraft shown. 

4.5 Wing Pressure Drag 

The calculations so far have defined those drag items which are inde- 
pendent of angle of attack and which compose the configuration's minimum drag 
level. The next step in determining the total drag is, to determine the 
variation of drag with lift and thereby define the polar shape. For the 

method herein described, polar shape is defined as a combination of the 

2 

theoretical induced drag (C / aR) and the corresponding value of wing 
pressure drag C^^^CQp) which varies as a function of wing lift and Mach 
number. 
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The calculation of wing pressure drag for the two representative 
configurations is begun by the tabulation of the required wing parameters 
from tables 2, 3, and 4: 



A-4F 

RA-5C 

^LdeS 

0.299 

0.365 

”des 

0.875 

0.889 

AR (t/c)^^^ 

1.23 

1.37 

(t/c)^^^ (1 +-i^) 

0.454 

0.368 


Using these values, the curves of figures 16 through 27 are Interpolated 
between for the appropriate value of AM and the resultant value of 
calculated as a function of AC as shown below. The corresponding value of 
wing lift coefficient, C^, is then determined by the expression: 




DES 


(19) 


These calculations are as follows: 


A-AF 



M = 0.80 

AM = 0.80 - 0.875 = -0.075 


<1 


^CDp 

(t/c)^^^ (1 +-^) 


-0.3 

t 

-0.001 

0.00050 

0.00023 

-0.2 

0.099 

0.00140 

0.00064 

-0.1 

0. 199 

0.00315 

0.00143 

-0.05 

0.249 

0.00498 

0.00226 

0 

0. 299 

0.00625 

0.00284 

0.05 

0.349 

0.01205 

0.00547 

0.1 

0.399 

0.02200 

0,00999 

0.2 

0.499 

0.05500 

0.02497 

0.3 

0. 599 

0.11000 

0.04994 
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RA-5C: 




Moj- 0.90AM - 0,90 - 

0.889 « 0.011 

- 2.0 AM - 2.0 - 0.889 - 1.111 

ACdp 

^Cdp 

ACdp 


(t/c)^^^ (1 + 

(t/c)^^^ (1 + 5^) 

-0.3 

0.065 

0.0005 

0.00018 

0.00258 

0.00095 

-0.2 

0.165 

0.00205 

0.00075 

0.0202 

0.00743 

-0.1 

0.265 

0.00465 

0.00171 

0.0465 

0.01711 

-0.05 

0.315 

0.00575 

0.00212 

0.0780 

0.02870 

0 

0.365 

0,01090 

0.00401 

0.1070 

0.03938 

0.05 

0.415 

0,0170 

0.00626 

0.1400 

0.05152 

0.1 

0.465 

0.0280 

0.01030 

0.1580 

0.05814 

0.2 

0.565 

0.0510 

0.01877 

- 


0.3 

0.665 

0.1020 

0.03754 

- 

- 


The computation of total drag-due-to-lif t is now accomplished by using 

the equation 

c_ -Ac + c//ttAR (20) 

P ^ 

A comparison of results from this procedure versus that obtained from the 
basic data for the A-4F and RA-5C are given in figures 61 and 62. Figures 63 
through 66 present similar comparisons for a representative sampling of the 
other study aircraft. 

4.6 Total Configuration Polar 


The total drag polar at the example Mach number is now computed for both 
configurations by the summation of the lift dependent and the minimum drag 
values as previously computed. This calculation is as follows: 

(Equations (2)^ (12)), 
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AA-4F 

RA - 5C 


- 0.80 




M - 0.90 

Mo, ■ 

2.0 


MIN 

- 0.01828 




\lN ■ 0-01787 

°HIN " 

0.03105 


Cj^^/itAR 


S 


2 

/rrAR 


Op 

"S 

s 

-0.001 

0.0 

0.00023 

0.01851 

0.065 

0.00036 

0.00018 

0.01841 

0.00095 

0.03236 

0.099 

0.00107 

0.00064 

0.01999 

0.165 

0.00232 

0,00075 

0.02094 

0.00743 

0.04080 

0.199 

0.00433 

0.00143 

0.02404 

0.265 

0.00599 

0.00171 

0.02557 

0.01711 

0.05415 

0.249 

0.00678 

0.00226 

0.02732 

0,315 

0.00847 

0.00212 

0.02846 

0.02870 

0.06822 

0.299 

0.00978 

0.00284 

0.03090 

0.365 

0.01137 

0.00«01 

0.03325 

0.03938 

0.08180 

0.349 

0.01332 

0.00547 

0.03707 

0.415 

0.01470 

0.00626 

0.03883 

0.05152 

0.09727 

0.399 

0.01741 

0.00999 

0.04568 

0.465 

0.01845 

0.01030 

0.04662 

0.05814 

0.10764 

0.499 

0.02724 

0.02497 

0.07049 

0.565 

0,02724 

0.01877 

0.06388 


- 

0.599 

0.03925 

0.04994 

0.10747 

0.665 

0.03774 

0.03754 

0.09315 

- 

- 


Comps.rl.son of thsss dsts to thst of the reference bsslc dsts is given by 
figures 67 and 68. 


4 • 7 BUFFET ONSET 


The ability to predict quickly the buffet onset characteristics of a 
configuration to within an acceptable tolerance level is of equal importance 
to the drag calculation herein described. To this end, the computer code 
(EDET) has been structured to Include the correlation of buffet onset infor- 
mation so that the lift coefficient for this occurance ) can be 

estimated. B.O. 

The calculation for buffet onset is accomplished in the following manner 
for the two representative aircraft. First, the following data must be 
assembled: 


Aircraft 

AR (1 + 

cos A ,, 
c/4 



^ES 

1 ^^ 

A-4F 

3.135 

0.8367 

0.1775 

0.299 

0.875 


RA-5C 

3.730 

0.7934 

0.1357 

0.365 


0.2127 
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Using figure 28 at the appropriate value of M and (t/c), the calculation 
proceeds as follows: 




A-4F 



RA 

-5C 


1 

AM 

\.o. ■ \es 

X 0 

h)ZS 

\.o. 

AH 

\)ES 

X © 

c - c 

\.o. \es 

%.o. 

0.40 

-0.475 

0.0680 

0.2548 

0.5538 

-0.489 

0.0695 

0.3268 

0.6918 

0.60 

-0.275 

0.0480 

0,1798 

0.4788 

-0,289 

0,050 

0.2351 

0.6001 



0.0279 

0.1045 

0.4035 

-0.089 

0,0340 

0.1598 

0.5248 

0.90 

0.025 

0.0158 

0.0592 

0.3582 

0.011 

0,0320 

0.1504 

0.5154 


These computed values are compared to those presented in the basic data of 
references 3 through 21 in figure 69. Also Included are similar data 
computed for a representative sampling of the other subject aircraft. 
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Figure 59. 
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Figure 61. - Drag-due-to-llf t comparison, A-4F. 
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Figure 63. - Drag-due-to-lif t comparison, T-2B and KA*3B. 
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Figure 65. - Drag-due- to-lift comparison, F 5A. 
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APPENDIX A 


FORM FACTOR GENERATION 

The geometry related form factor for the fuselage, as presented in 
figure 6 of Section 1, was generated using results of a NASA test (Refer- 
ence 39). These data, which are representative of transport bodies, were 
used as the base for determining the minimum drag of a series of bodies of 
varying fineness ratio. These test results are compared on figure 70 with 
the estimated skin friction drag. A decided increase in pressure drag with 
reduced body fineness ratio is evidenced. On figure 71 these same data are 
presented as a form factor and compared with fuselages defined as an equiva- 
lent length to diameter ratio. Included on this figure are the fuselage form 
factors that would be predicted from other sources. 

The wing form factors for advanced and conventional sections, as given 
by figure 7, are based on information contained in references 34 and 40. 

On figure 72, this form factor has been developed for advanced airfoil sec- 
tions and data for a NACA 65 series section is also presented for reference. 
The parameter /^Cp ■ FF is the 2-dimensional minimam section drag co- 
efficient divided by twice the flat plate skin friction coefficient at the 
test Reynolds number. The factor of 2 accounts for both upper and lower 
surfaces. These data are from 2-D tests on a NACA 65, 213 a » 0.5 airfoil, 

9 percent state-of-the-art airfoil, 10 percent, 11 percent, 12 percent, and 
21 percent thickness advanced airfoils respectively. Flagged versus 
unflagged symbols represent the same model tested in two different facilities. 
On figure 73, the average fairing of C /2CF is noted as the form factor 
versus section thickness ratio. The conventional line, which was derived 
from RAS Dats Sheets (reference 41) assuming transition at the leading edge, 
is confirmed by results from the NACA 65 series airfoil. At a thickness 
ratio of 10 percent the advanced airfoil appears to carry an approximate 

10 percent increased subsonic pressure drag over the conventional airfoil 
sections. 


160 


ORIGINAL PAGE IS 
OF. POOR QUALITY 



X61 

































(AIRF01D/2 Cc (FLAT PLATE) 


OR/GINAL page is 
Pf POOR 






Figure 72. - Ratio of minimum drag to theoretical skin friction drag for 

conventional, state of the art, and advanced airfoil sections. 
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Figure 73. - Wing section form factor correlation. 
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